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Resumen 


Paredes-Trejo, F., Guevara-Pérez, E., Barbosa-Alves, H., 
Uzcátegui-Briceño, C. (noviembre-diciembre, 2015). Tendencia 
de la precipitación estacional e influencia de El Niño — Oscilación 
Austral sobre la ocurrencia de extremos pluviométricos en la 
cuenca del lago de Valencia-Venezuela. Tecnología y Ciencias del 
Agua, 6(6), 33-48. 


La cuenca del lago de Valencia (CELV) es la cuenca endorreica de 
mayor tamaño en Venezuela. Por su elevada densidad poblacional 
e industrial es susceptible a los extremos pluviométricos. Se sabe 
que el fenómeno ENOA (El Niño-Oscilación Austral) modula 
las lluvias en el territorio venezolano, pero no se ha explorado 
su incidencia en detalle en la CELV. En este estudio se analiza 
la tendencia espacial y temporal de la precipitación estacional 
y se explora la asociación entre la ocurrencia de meses con 
extremos pluviométricos y las fases de ENOA (El Niño, La Niña, 
neutro) en la CELV. Se seleccionaron ocho estaciones climáticas 
con buena calidad de registros. Los periodos 1934-2005 y 1966- 
1992 se adoptan para los análisis a escalas local y regional. Se 
identificaron los meses de la temporada seca y húmeda. En 
cada estación se calculó la precipitación acumulada estacional 
y se evaluó su tendencia de largo plazo utilizando la prueba 
de Mann-Kendall. Se categorizó la precipitación mensual local 
y estacional en extrema seca (ES), no extrema (NE) y extrema 
húmeda (EH), usando como umbrales los percentiles 10 y 90. Se 
analizó la ocurrencia probabilística espacial y simultaneidad de 
un mes ES, NE y EH, según la temporada y fases de ENOA. La 
asociación entre ENOA y la precipitación estacional se explora 
con una prueba Chi-Cuadrado. Se encontró lo siguiente: no 
existen tendencias locales de largo plazo en la precipitación total 
estacional; la ocurrencia de extremos pluviométricos estacionales 
está parcialmente asociada con los eventos El Niño/La Niña; la 
incidencia de extremos pluviométricos podría estar vinculada 


con factores climáticos locales. 
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Abstract 


Paredes-Trejo, F., Guevara-Pérez, E., Barbosa-Alves, H., € Uzcátegui- 
Briceño, C. (November-December, 2015). Seasonal Trend of the Rainfall 
and the Influence of El Niño-Southern Oscillation on the Occurrence 
of Extreme Rainfalls at the Watershed of Valencia's Lake, Venezuela. 
Water Technology and Sciences (in Spanish), 6(6), 33-48. 


The Valencia lake basin is the endorheic watershed larger of Venezuela 
(VLB). VLB is densely populated and industrialized; therefore the 
extreme climate events can cause severe impacts. In the Venezuelan 
territory, it is known that ENSO (El Niño-Southern Oscillation) can 
modulate the rainfalls, however no has been explored in detail as ENSO 
can affect the rainfalls on VLB. This study analyzes the spatio-temporal 
trends of the seasonal rainfalls and explores the association between 
the occurrence of extreme rainfall months and the phases of ENSO (El 
Niño, La Niña and Neutral) on VLB. Eight stations were selected by 
the quality control of their registers. We considered two periods for our 
analysis: 1934-2005 (local scale) and 1966-1992 (regional-local scale). 
In addition, were identified the months of occurrence of the rainy and 
dry seasons. The accumulated rainfall by season in each station was 
calculated and after was explored the occurrence of a long-term trend 
by Mann-Kendall test. Seasonal precipitation in each station and season 
was categorized in extreme dry, extreme-no, or extreme wet (ED, EN 
and EW) using the 10th and 90th percentiles as threshold. The likelihood 
of occurrence of a month ED, EN or EW at local scale according to 
season and phase of ENSO was estimated. The association between these 
categorical variables was analyzed by chi-square test for independence. 
The results more remarkable were: there is no evidence of a long-term 
trend at local scale on seasonal precipitation; the Niño/Niña episodes are 
partially associated with the occurrence of extreme rainfall at seasonal 
scale; a high-proportion of the extreme rainfall events could have been 
driven by local factors which were not evaluated here. 


Keywords: Valencia's Lake, El Niño-Southern Oscillation, climate 
change, extreme rainfall. 
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Introducción 


En el Cuarto Informe de Evaluación del Grupo 
Intergubernamental de Expertos sobre el Cambio 
Climático se indicó que entre los años 1900 y 
2000, la precipitación anual aumentó de forma 
notable en la región ubicada al norte de América 
del Sur (Grupo Intergubernamental de Expertos 
sobre el Cambio Climático, 2007). En esa región, 
el fenómeno meteorológico El Niño-Oscilación 
Austral o ENSO (ENOA) es uno de los com- 
ponentes del sistema climático global que más 
influye sobre la variabilidad pluviométrica (Ar- 
ntz £ Fahrbach, 1991; Giddings $ Soto, 2006). 
ENOA ha sido ampliamente relacionado con las 
impactos climáticos más importantes ocurridos 
en Centroamérica, México y el Cono Sur (Banco 
Mundial-Naciones Unidas, 2001; Comisión Eu- 
ropea, 2009). Esta influencia se extiende a los 
países situados dentro y en la periferia de las 
cuencas del Pacífico e Índico (Allan, Lindesay, 
£z Parker, 2000), el Trópico y ciertas regiones 
extra-tropicales de Norteamérica (Magaña éz 
Ambrizzi, 2005; PNUMA-CEPAL, 2010). 
Venezuela es un país tropical situado en el 
extremo norte de Sudamérica (- 1-12 N, 60-742 
O), con la precipitación sobre este territorio 
principalmente modulada por la actividad y 
ubicación de la Zona de Convergencia Intertro- 
pical (ZCIT), en la cual confluyen los vientos 
alisios del noreste y del sudeste, estructurando 
un subsistema de baja presión que se desplaza 
en sentido meridional (Goldbrunnex, 1984; 
Pulwarty, Barry, Hurst, Sellinger, £ Mogollon, 
1998). La ZCIT es parte del sistema climático 
global; por tanto, responde rápidamente a las 
variaciones en la circulación atmosférica tropical 
(Martelo, 2003a). En el contexto espacial, la pre- 
cipitación es heterogénea debido al control que 
ejercen las barreras orográficas sobre los vientos 
de superficie (Insel, Poulsen, 4 Ehlers, 2010). 
Mientras que en el contexto temporal, varios 
modos de variación han sido identificados en las 
series de precipitación disponibles (Winemiller, 
1990). A escala decadal, por ejemplo, Martelo 
(2003b) notó que los promedios pluviométricos 
revelaban la ocurrencia de un patrón débil- 
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mente alternante; así, en la década 1951-1960 
prevaleció una condición normal; 1961-1970 fue 
lluviosa; 1971-1980, seca; 1981-1990, lluviosa, y 
la década incompleta 1991-1998, seca. 

La variabilidad pluviométrica sobre el terri- 
torio venezolano en el contexto temporal puede 
ser modulada por las fases de ENOA tanto cáli- 
da o El Niño como fría o La Niña (Pulwarty, Ba- 
rry, € Riehl, 1992; Cárdenas, García, €: Gil, 2002; 
Paredes, Guevara, Uzcategui, $: Garbi, 2008). El 
efecto de ENOA sobre el territorio venezolano 
no es homogéneo (C?NMeH-CONICIT, 1998); 
así, en las regiones de Guayana y occidente, El 
Niño tiende a propiciar la ocurrencia de veranos 
más secos que lo normal y un retardo en el inicio 
de la temporada lluviosa; mientras que en las re- 
giones central, llanera y oriental, la influencia es 
mucho menos evidente, encontrándose casos en 
los que los años El Niño coinciden con inicios de 
temporada de lluvia temprana, normal o tardía. 
En contraste, un ENOA neutro, es decir, océano 
Pacífico con condiciones térmicas superficiales 
cercanas al promedio histórico, suele coincidir 
con un inicio temprano de la temporada lluvio- 
sa. Cabe destacar que en Guayana, La Niña se 
asocia con temporadas secas más lluviosas de 
lo normal y por tanto la ocurrencia de caudales 
relativamente altos en el río Caroní, mientras 
que en los Llanos se vincula con una entrada 
tardía de la temporada de lluvias (Rogers, 1998; 
Martelo, 2003b). 

A gran escala, la influencia de ENOA sobre 
el régimen de las precipitaciones en el territorio 
venezolano ha sido ampliamente estudiada por 
Caviedes (1998), Cárdenas y De Grazy (2003), 
Ministerio del Ambiente y Recursos Naturales 
Renovables (2007), y Rollenbeck y Anhuf (2007); 
sin embargo, las principales características 
espaciales y temporales de esta modulación 
sólo se han analizado en suficiente detalle en 
algunas regiones del país (Guevara €: Paredes, 
2007; Pierre €: Tirado, 2007; Pérez, 2012). Una 
de las regiones donde no se conoce muy bien la 
incidencia de ENOA sobre la precipitación es la 
cuenca del lago de Valencia. El relieve dominan- 
te en esta cuenca incrementa la vulnerabilidad 
ante la ocurrencia de extremos pluviométricos. 
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En lo sucesivo, el término extremo pluviométri- 
co se refiere a un episodio anómalamente seco o 
húmedo, con una persistencia de al menos dos 
meses consecutivos. 

El lago de Valencia es el cuerpo de agua dul- 
ce endorreico más grande de Venezuela (Díaz, 
2006), ocupa una depresión tectónica denomi- 
nada Graben de Valencia entre la Cordillera de 
la Costa al norte y la Serranía del Interior al sur 
(Bradbury et al., 1981). Está rodeado por una 
región altamente industrializada y densamente 
poblada (Sequera, 1994). A grandes rasgos, los 
suelos son destinados a actividades agrícola 
(34%), urbana (18%) e industrial (3%). El cuerpo 
de agua del lago tiene 359 km” de superficie, 
6.30 km* de volumen y 18 m de profundidad, al 
cual tributan 18 subcuencas (Guevara €: Már- 
quez, 2012; Dourojeanni, Jouravlev, $: Chávez, 
2002). Los ríos con régimen permanente son 
Guey, El Limón, Las Delicias, Turmero, Aragua, 
Cabriales, Los Guayos, Nepe y Guacara, junto al 
Caño Central y los canales Papelera, Sudantex 
y Corpoindustria (Filippone, 1999). El escurri- 
miento superficial medio hacia el lago es de 9 
nY/s durante la temporada seca y de 19 n?*/s 
durante el periodo húmedo (Guevara, Guevara, 
éz García, 2008). 

Durante la última década, el nivel del lago 
de Valencia ha aumentado en forma sostenida 
debido al transvase desde los embalses Pao Ca- 
chinche y Pao La Balsa, que alimenta el Sistema 
Regional del Centro 1 (principal acueducto del 
eje Valencia-Maracay); el desvío del río Cabria- 
les en el año 1979 hacia el lago, y el vertido con- 
tinuo de aguas servidas urbanas e industriales. 
Entre 1995 y 2000, el nivel de ascenso del espejo 
de agua del lago subió a una velocidad media 
de 0.40 m/año, alcanzando un máximo relativo 
de 413.36 msnm en el año 2012 (Ministerio del 
Ambiente y Recursos Naturales Renovables, 
2007). Durante la temporada de lluvias del año 
2013, la Comisión de Saneamiento y Control 
de Nivel del Lago de Valencia (dependencia 
gubernamental) reportó la ocurrencia de una 
cota cercana a los 414 msnm y varios impactos 
de importancia: inundación parcial de algunos 
sectores situados al sur del municipio Girardot 


del estado Aragua, como Mata Redonda y La 
Punta; interrupción temporal del acceso vial a la 
Isla La Culebra, y anegamiento permanente de 
las bases del tramo vial Valencia-Gúigúe, entre 
otros. 

La vocación agrícola de la cuenca del lago 
de Valencia, especialmente en las tierras bajas, 
la hace muy vulnerable a los episodios secos 
persistentes. Por ejemplo, algunos cultivos de 
secano como sorgo industrial, hortalizas, frutas 
y cítricos merman sus rendimientos durante la 
ausencia prolongada de lluvias, lo que afecta al 
sector agrícola de la región (Marín, 2002). 

En las referencias se encuentran estudios 
que caracterizan la distribución espacial de 
las lluvias en la cuenca del lago de Valencia 
(Guevara et al., 2008); otros que estiman algunos 
parámetros hidrológicos relevantes (Ascenzi, 
Mora, €: Pino, 2007) y que evalúan bajo un 
enfoque probabilístico la precipitación anual 
y estacional (Paredes et al., 2013). Sin embargo, 
se desconoce si la precipitación en la cuenca ha 
aumentado o disminuido en el tiempo; tampoco 
se sabe si las temporadas anómalamente secas O 
húmedas están asociadas con las fases de ENOA 
o si los extremos pluviométricos tienden a ser 
más frecuentes. Por ello, el propósito de esta 
investigación es la de explorar la influencia 
del fenómeno meteorológico ENOA sobre los 
extremos pluviométricos en la mencionada 
cuenca. En el estudio se plantean dos objeti- 
vos relacionados con el análisis de las series 
pluviométricas mensuales disponibles en un 
conjunto de estaciones situadas dentro de la 
cuenca: a) analizar la tendencia de largo plazo 
del acumulado pluviométrico estacional, y b) 
analizar la vinculación entre las fases de ENOA 
y los extremos pluviométricos estacionales bajo 
un enfoque probabilístico. 


Metodología 
Área de estudio 
La cuenca del lago de Valencia se ubica en la 


Región Central de Venezuela. Cubre una super- 
ficie de 2 943 kmY, distribuida entre los estados 
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Carabobo (53.13%) y Aragua (46.87%). Se extien- 
de entre las latitudes norte 9” 55' 4.26” y 10” 24” 
41.74”, longitudes oeste -68” 4” 21.2” a -67* 16' 
30.15”, como se aprecia en la figura 1. Tiene un 
rango hipsométrico que va de 426 msnm (nivel 
medio del lago) a 2 439 msnm (Pico Codazzi). 


Registros pluviométricos 


Se usaron las series pluviométricas mensuales 
de las estaciones climáticas localizadas en el 
área de estudio y administradas por el Instituto 
Nacional de Meteorología e Hidrología (http: / / 
www.inameh.gob.ve/). Se seleccionaron las 
series que cuentan con una longitud de registro 
mensual (LRM) igual o mayor que 30 años. 
Con base en este criterio se eligieron la serie 
de registros de ocho estaciones pluviométricas 
con 31 años > LRM > 59 años, con las siguientes 
características: registros discontinuos, extensión 
media de 42.50 años y mediana de la extensión 


de 42.50 años. Estos periodos se enmarcan en 
una ventana temporal discontinua que va desde 
el año 1934 hasta el año 2005. En lo sucesivo, a 
este grupo se denomina muestra. 


Fases del procedimiento 


1. Evaluación de la calidad de las series. En la 
muestra, la estacionariedad se evaluó utili- 
zando el estadístico t (criterio de aceptación, 
p < 0.05), tomando como variable de análisis 
la pendiente de la recta que relaciona la 
precipitación anual y el tiempo en años; la 
autocorrelación serial se evaluó con la prue- 
ba de Durbin-Watson (Fox, 2008); los años 
no homogéneos se detectaron con la prueba 
de Easterling-Peterson (Buishand, 1982). Las 
ocho series de registros pluviométricos usa- 
dos pasaron las pruebas mencionadas y, por 
tanto, se usaron en los análisis posteriores 


(figura 1). 
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Figura 1. Cuenca del lago de Valencia. Nota: 1, río Cabriales; 2, río Los Guayos; 3, río Guacara; 4, río Mariara; 5, río Tapatapa; 


6, río Turmero; 7, laguna de Taguayguay; 8, embalse Zuata; 9, río Aragua; 10, río Las Minas; 11, río Tocorón; 12, río Gúigiey; 


13, Caño Central; 14, río Las Dos Bocas; 15, quebrada Honda; 16, río Maracay; 17, Dique Las Delicias. Los triángulos señalan la 


ubicación de las ocho estaciones usadas; el número es el serial de la estación, según la nomenclatura del Instituto Nacional de 


Meteorología e Hidrología venezolano. 
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2) Análisis de la tendencia anual en los montos 
pluviométricos estacionales. En principio 
se identificaron los meses de inicio y fi- 
nalización de las temporadas húmeda y 
seca, tomando como criterio el signo del 
coeficiente pluviométrico a escala mensual 
(Carillo, 1999), resultando: noviembre-abril, 
temporada seca; mayo-octubre, temporada 
húmeda. En un año se consideró que la 
temporada seca va de enero hasta abril para 
evitar el solapamiento de años consecutivos. 
Se determinó en cada serie y año la precipi- 
tación de la temporada (PT), el coeficiente 
de variación de la PT, monto y año de ocu- 
rrencia de la PT máxima y mínima. Bajo este 
enfoque, PT es la lluvia acumulada anual 
durante la temporada analizada (seca o 
húmeda). Por último, se aplicó la prueba de 
Mann-Kendall basada en rangos (Kendall, 
1975), para identificar tendencias de largo 
plazo estadísticamente significativas (p < 
0.05) en la PT. Las variables PT durante las 
temporadas seca y húmeda fueron mapea- 
das para identificar la ocurrencia de algún 
patrón espacial. 


3) Análisis temporal de las anomalías pluvio- 
métricas extremas estacionales. Se calculó 
el periodo común de registros entre las 
series que conforman la muestra, resultan- 
do 1966-1992. En cada serie se transformó 
la precipitación mensual (P,) en anomalía 
mensual (a); para ello se calcularon los 
percentiles 10 (P,,) y 90 (P,,) de cada mes 
i (i=1... 12), considerando todos los años ¡ 
del periodo común (¡= 1966... 1992), luego 
se codificó cada observación como sigue: si 
A (mes extremo seco o ES); si 
PE, >4,=*+1 (mes extremo húmedo o 
EA) siP,< PS Po => a,= 0 (mes no extremo 
o NE). Se evaluó la tendencia estacional de 
los meses extremos (ES, EH), aplicando la 
prueba de Mann-Kendall al vector resul- 
tante de parear, por separado, el número de 
meses ES y EH con los años disponibles. Las 
tendencias estacionales (pendientes) fueron 
mapeadas para identificar la ocurrencia de 
algún patrón espacial. 


4) Relación probabilística entre los niveles plu- 
viométricos estacionales en cada estación y 
las fases de ENOA. Se calculó la probabili- 
dad mensual de ocurrencia de un mes ES, 
EH y NE en cada estación, temporada y fase 
ENOA. La fase de ENOA que prevaleció en 
un mes determinado se tomó de la categori- 
zación de ENOA disponible en el portal web 
del Centro de Predicción Climática de la 
Agencia Nacional del Océano y Atmósfera 
de los Estados Unidos (www.cpc.ncep.noaa. 
gov). Dado que la precipitación estacional 
en cada estación y las fases ENOA son 
analizadas como factores, se estructuró una 
tabla de contingencia de dos vías, sobre la 
cual se aplicó una prueba x? para evaluar la 
asociación entre ambas variables categóricas 
a un nivel de confianza de 95% (Greenwood 
éz Nikulin, 1996). Para estimar de modo in- 
directo la cobertura espacial promedio de 
los meses ES, EH y NE en cada temporada 
durante el periodo 1966-1992, se calculó 
el promedio porcentual de estaciones que 
mostraron un mes ES, EH y NE en cada año 
y temporada; este resultado fue agrupado 
por fase ENOA observada, para evaluar la 
asociación entre estas variables. 


Herramientas computacionales utilizadas en 
los análisis 


La calidad de las series se evaluó con el software 
Anclim (Stépánek, 2008). La conversión de las 
series mensuales pluviométricas a anomalías, 
el análisis probabilístico y estadístico fueron 
realizados en ambiente R 7. 3.2.1. El mapeado 
se hizo con SAGA-GIS 2.1.4'M, 


Resultados 


Tendencia de los montos pluviométricos 
estacionales en el periodo no común 1934-2005 


Las series pluviométricas analizadas tienen 
diferentes fechas de inicio y fin dentro de la 
ventana temporal 1934-2005. El análisis de cada 
serie por separado, considerado todos sus regis- 
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tros, arrojó lo siguiente: la precipitación media 0.026 eventos ES /año; en el grupo de EH, -0.009 
regional durante la temporada seca (enero-abril) eventos EH /año. Entre mayo y octubre (tempo- 
es igual a 86.75 mm (cuadro 1); el mínimo y rada húmeda), la tendencia de mayor magnitud 
máximo local ocurren en Maracay (50.17 mm) en el grupo de ES se registró en la Colonia El 
y Las Dos Bocas (141.75 mm), respectivamente. Trompillo (-0.032 eventos ES / año); en el grupo 
El coeficiente de variación regional promedio es de EH ocurrió en Guacara (-0.038 eventos EH / 
76.14%, variando de 62.22% en Las Dos Bocas año). Globalizando, el grupo de ES muestra, en 
a 86.95% en Maracay. Durante la temporada promedio, -0.003 eventos ES /año; mientras que 
húmeda (mayo-octubre), el promedio regional el grupo de EH, -0.006 eventos EH /año. Ningu- 
de precipitación es igual a 873.12 mm, con na de estas tendencias resultó estadísticamente 
un mínimo y máximo local, ubicados en Las significativa (p < 0.05). 
Cenizas (715.67 mm) y Las Dos Bocas (1234.87 
mm), respectivamente (figura 1); el coeficiente Característica espacial y temporal de los 
de variación regional promedio es de 20.13%; extremos pluviométricos estacionales según la 
las tendencias de la precipitación en las tempo- fase de ENOA prevaleciente durante el periodo 
radas húmeda y seca (expresadas como mm/ 1966-1992 
año) fueron estadísticamente no significativas 
en todas las estaciones (p < 0.05). El cuadro 3 muestra la probabilidad ocurrencia 
de un mes ES, EH o NE durante las temporadas 
Tendencia de las anomalías pluviométricas seca y húmeda (mayo-octubre) para cada esta- 
extremas estacionales en el periodo común ción dentro de la cuenca del lago de Valencia 
1966-1992 en el periodo 1966-1992. A manera de ejemplo, 
cuando está presente un episodio El Niño a 
El cuadro 2 muestra la tendencia anual de ocu- lo largo de la temporada seca (enero-abril), 
rrencia de meses extremos en las temporadas la estación 417 registró una probabilidad de 
secas y húmedas para el periodo 1966-1992. ocurrencia de un mes ES, NE y EH igual que 
En general, durante la temporada seca (enero- 31, 56 y 13%, respectivamente. Ante un episodio 
abril), la tendencia de mayor magnitud en el La Niña, estas probabilidades son de 13, 71 y 
grupo de meses extremos secos (ES) ocurre 17%, respectivamente. Mientras que en la fase 
en Las Cenizas (0.046 eventos ES/año); en el neutra de ENOA se dan valores de 28, 63 y 9%, 
grupo de meses extremos húmedos (EH) se respectivamente. La asociación entre ENOA y 
observa en Santa Cruz (-0.020 eventos EH/año) los niveles pluviométricos (ES, NE y EH) en 
y Guacara (-0.020 eventos EH/año). En cuanto la temporada seca resulta estadísticamente 
a la temporada seca, el grupo de ES promedia significativa sólo en las estaciones 489 y 491 (p 


Cuadro 1. Montos pluviométricos estacionales. Periodo no común 1934-2005. 


Serial ión Periodo Estación seca Estación húmeda 
(años) (P(mm)/CV(%)) (P(mm)/CV(%)) 

417 Santa Cruz 1966-1999 70.36/72.95 840.09 /19.81 
452 Guacara 1949-1993 70.66/78.30 782.18/19.51 
466 Maracay 1934-1992 50.17/86.95 774.42 /22.12 
488 Colonia El Trompillo 1960-2005 79.52/86.09 858.49 /24.42 
489 Agua Blanca 1934-2005 119.69 /73.74 1 014.78/22.52 
491 Las Dos Bocas 1949-2005 141.75/62.22 1 231.88/13.38 
497 Las Cenizas 1960-2003 91.19/73.59 715.67 /22.71 
1494 Embalse Taiguaiguay 1951-1999 70.63/75.31 766.60 /16.58 
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Cuadro 2. Tendencia anual de ocurrencia de meses extremos en las temporadas secas y húmedas. Periodo común 1966-1992. 


Serial on Temporada seca Temporada húmeda 
(ES/EH) (eventos/año) (ES/EH) (eventos/año) 
417 Santa Cruz 0.0031 /-0.0195 0.0024 /-0.0122 
452 Guacara 0.0250/-0.0195 -0.0171/0.0379 
466 Maracay 0.0269 /-0.0006 -0.0073/0.0214 
488 Colonia El Trompillo 0.0140/-0.0183 -0.0324/0.0024 
489 Agua Blanca 0.0379 /0.0018 0.0098 /0.0244 
491 Las Dos Bocas -0.0037 /0.0122 0.0116/-0.0018 
497 Las Cenizas 0.0458 / -0.0073 0.0037 /-0.0092 
1494 Embalse Taiguaiguay 0.0317 /-0.0183 0.0031 /-0.0183 


Nota: ES, mes extremo seco; EH, mes extremo húmedo. 


< 0.05), ambas localizadas en la zona media-alta 
de la cuenca y en el flanco sur del lago (figura 2). 
Durante la temporada húmeda (mayo-octubre), 
sólo en la estación 1494, la prueba de indepen- 
dencia basada en el estadístico x? sugiere una 
asociación entre la actividad de ENOA y la 
precipitación estacional (cuadro 4). 

El cuadro 3 muestra la asociación cuantita- 
tiva entre las fases de ENOA (neutro, El Niño, 
La Niña) y la señal pluviométrica dominante 
en cada estación y temporada (ES, EH o NE) 
bajo un enfoque probabilístico, pero no brinda 
información sobre la extensión espacial de estos 
eventos a escala estacional. La figura 3 com- 
pensa esta limitación debido a que muestra el 
promedio regional porcentual de estaciones en 
los que se observó un mes ES, EH o NE durante 
las temporadas seca (enero-abril) y húmeda 
(mayo-octubre) a lo largo del periodo 1966-1992, 
junto a la fase ENOA registrada durante cada 
temporada (neutro, El Niño, La Niña). 

Un promedio porcentual alto en el nivel ES, 
EH o NE indica que esta señal predominó en 
el grupo de estaciones. Por ejemplo, la figura 
3 revela que en la temporada húmeda del año 
1973 (panel b), cerca de 70% de las estaciones 
tienen un mes ES junto a una condición La 
Niña (condición seca generalizada); mientras 
que en la temporada seca del año 1981 (panel 
a), alrededor de 63% de la estaciones muestra un 
mes EH acoplado a una condición ENOA neutra 
(condición húmeda generalizada). 


Discusión 


Tendencia de los montos pluviométricos esta- 
cionales en el periodo analizado 


Desde el punto de vista espacial, los montos 
pluviométricos estacionales tienden a dismi- 
nuir conforme aumenta la cercanía al lago de 
Valencia y el terreno se hace plano (cuadro 1; 
figura 2). Los valores más bajos se observan 
en las planicies situadas al sur de la Troncal 5 
(tramo Los Guayos-Maracay) y al norte de la 
intercomunal Flor Amarilla-Magdaleno (flanco 
sur del lago). Los más altos se dan en el extremo 
suroeste (Cerro Las Dos Bocas). 

Los registros eólicos de las estación Nagua- 
nagua (10? 14 58” N y 68 1” 0” O), situada en 
un extenso valle bordeado por las Filas, La Gua- 
camaya y Orégano (ambas alineadas en sentido 
norte-sur), muestra una interesante característi- 
ca: los vientos de superficie durante la tempora- 
da seca provienen, en orden de ocurrencia, del 
NO, NNO y ONO, con una velocidad promedio 
de 16.2 km/h. En la temporada húmeda amino- 
ran su velocidad hasta unos 10 km/h y cambian 
su dirección hacia el NO, ENE y SO (en orden 
de ocurrencia). En este punto, cabe destacar que 
sobre el territorio venezolano, el flujo eólico su- 
perficial a gran escala y a lo largo del año tiene 
direcciones prevalecientes E, ENE y NE debido 
a la influencia de los alisios sobre el norte de Su- 
damérica (González-Longatt, 2015). Sin embar- 
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Cuadro 3. Probabilidad de ocurrencia de un mes extremo seco (ES), no extremo (NE) o extremo húmedo (EH) durante las 


temporadas seca y húmeda en las estaciones de la cuenca del lago de Valencia de acuerdo con la fase ENOA observada (neutro, 


El Niño, La Niña). Periodo común 1966-1992. La probabilidad de ocurrencia se expresa como una fracción unitaria. 
Temporada seca (enero-abril) Temporada húmeda (mayo-octubre) 
Serial: 417 Serial: 417 
Condición ES NE EH ES NE EH 
Neutro 0.28 0.63 0.09 0.11 0.77 0.11 
El Niño 0.31 0.56 0.13 0.08 0.75 0.17 
La Niña 0.13 0.71 0.17 0.17 0.79 0.04 
Serial: 452 Serial: 452 
Condición ES NE EH ES NE EH 
Neutro 0.33 0.53 0.14 0.10 0.80 0.09 
El Niño 0.25 0.65 0.10 0.10 0.80 0.10 
La Niña 0.42 0.54 0.04 0.14 0.69 0.17 
Serial: 466 Serial: 466 
Condición ES NE EH ES NE EH 
Neutro 0.48 0.39 0.13 0.09 0.78 0.13 
El Niño 0.35 0.55 0.10 0.23 0.70 0.07 
La Niña 0.25 0.67 0.08 0.06 0.83 0.11 
Serial: 488 Serial: 488 
Condición ES NE EH ES NE EH 
Neutro 0.20 0.67 0.13 0,13 0.78 0.09 
El Niño 0.25 0.65 0.10 0.10 0.80 0.10 
La Niña 0.21 0.71 0.08 0.08 0.75 0.17 
Serial: 489* Serial: 489 
Condición ES NE EH ES NE EH 
Neutro 0.13 0.73 0.14 0.11 0.78 0.10 
El Niño 0.35 0.50 0.15 0.07 0.80 0.13 
G La Niña 0.13 0.88 0.00 0.14 0.75 0.11 
de Serial: 491* Serial: 491 
5 Condición ES NE EH ES NE EH 
3 Neutro 0.20 0.63 0.17 0.11 0.77 0.11 
hs El Niño 0.05 0.95 0.00 0.10 0.77 0.13 
FE La Niña 0,17 0.79 0.04 0.11 0.81 0.08 
3 Serial: 497 Serial: 497 
5 Condición ES NE EH ES NE EH 
E Neutro 0.16 0.70 0.14 0.07 0.80 0.13 
a El Niño 0.40 0.50 0.10 0.20 0.77 0.03 
E La Niña 0.25 0.71 0.04 0.14 0.72 0.14 
= Serial: 1494 Serial: 1494 * 
ES Condición ES NE EH ES NE EH 
S Neutro 0.27 0.63 0.11 0.07 0.84 0.08 
a El Niño 0.35 0.55 0.10 0.20 0.77 0.03 
z La Niña 0.25 0.63 0.13 0.14 0.61 0.25 
3 Nota: *la ocurrencia de un mes extremo seco, no extremo o extremo húmedo en la estación indicada no es independiente de las condiciones 
3, ENOA, con un nivel de confianza de 95%. 
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Figura 2. Modelo digital de elevación de la cuenca del lago de Valencia. Nota: los triángulos señalan la ubicación de las ocho 


estaciones usadas; el número es el serial de la estación, según la nomenclatura del Instituto Nacional de Meteorología e 


Hidrología. 


go, estas direcciones suelen ser poco frecuentes 
en los valles intermontanos profundos debido 
a que el flujo eólico es reorientado (Golbrunner, 
1984; Foghin-Pillin, 2002). La diferencia entre 
los patrones eólicos dominantes durante las 
temporadas seca y húmeda, observada en la 
estación Naguanagua, sugiere que las barreras 
orográficas encauzan las corrientes de aire hacia 
el fondo de los valles, donde interactúan con las 
circulaciones eólicas locales (sistema lago-brisas 
del valle-montaña), favoreciendo o inhibiendo 
la formación de lluvias, según la dirección del 
flujo dominante y su contenido de humedad. El 
control que ejerce la orografía sobre la distribu- 
ción de las lluvias en las cuencas de montañas 
ha sido bien documentado (Beniston, Diaz, € 
Bradley, 1997; Guan, Huang-Hsiung, Makhnin, 
Xied, £ Wilsone, 2009; Brito et al., 2010) 
Dentro de la cuenca del lago de Valencia se 
distinguen dos núcleos húmedos que permane- 
cen en esa condición durante las temporadas 
seca y húmeda. Uno de ellos se localiza en el 


río Las Dos Bocas, entre las Filas, Los Aguaca- 
tes y Toma Las Manos; una región montañosa 
ubicada al suroeste de la cuenca (figura 2). El 
segundo núcleo se sitúa al noreste de la línea 
divisoria entre las cuencas del lago de Valencia 
y el río Chuao (lugar denominado Topo El Gua- 
yabo). El primero de estos núcleos fue notado 
previamente por Paredes et al. (2013), quienes 
advirtieron el carácter permanente de su plu- 
viosidad a lo largo del año. Estos núcleos tienen 
en común el hecho de estar al fondo de valles 
estrechos, flanqueados por altas montañas, y 
próximos al parteaguas o divisoria. La figura 
2 revela que la cuenca presenta una amplia red 
de profundos valles al norte y sur del lago, por 
donde el viento de superficie podría ser canali- 
zado desde el lago (subsistema climático local) o 
cuencas hidrográfica vecinas hacia la zona más 
elevada. En esta circunstancia, las corrientes de 
aire ascenderían y enfriarían progresivamente 
hasta alcanzar el nivel de condensación sobre 
la vertiente ascendente (ladera Barlovento), fa- 
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Figura 3. Promedio regional de ocurrencia de un mes extremo seco (ES), no extremo (NE) y extremo húmedo (EH) durante 


las temporadas de seca (panel a) y húmeda (panel b) en la cuenca del lago de Valencia en el periodo 1966-1992. En el eje 


superior de cada panel se indica la fase ENOA observada (neutro, El Niño o La Niña). Nota: el año 1992 fue omitido porque las 


estaciones tienen datos faltantes en las temporadas seca y húmeda. 


voreciendo la formación de nubes generadoras 
de lluvias. Este hipotético mecanismo permite 
explicar los altos valores de precipitación es- 
tacional y su persistencia en Las Dos Bocas y 
Topo El Guayabo. Bajo este enfoque, el viento 
superficial que llega a Las Dos Bocas es probable 
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que ingrese siguiendo la ruta Pirital-Boquerón- 
El Yagual (corredor intramontano). 

En el contexto temporal, la evidencia obser- 
vacional para los periodos 1934-2005 y 1966- 
1992 no indica la ocurrencia de una disminución 
o aumento sostenido en los totales pluviométri- 
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cos estacionales extremos (cuadro 2). En general, 
esto sugiere que las variaciones pluviométricas 
estacionales en la cuenca del lago de Valencia en 
estos periodos fueron temporalmente estacio- 
narias, aunque es evidente que la distribución 
espacial de las lluvias es controlada en gran 
medida por la orografía. 

El contraste de la figura 1, y los cuadros 1 y 
2 revela que durante la temporada seca se da 
un ligero incremento en la ocurrencia de meses 
ES en la cuenca alta del río Las Minas (tribu- 
tario de la laguna de Taguaguay), y una leve 
disminución de meses EH en las nacientes del 
río Guacara y en el suroeste de Cagua (llanuras 
del río Aragua). Mientras que en la temporada 
húmeda se nota una sutil disminución en la 
ocurrencia de meses ES en la cuenca media 
del río Gúiigite y meses EH en las nacientes del 
río Guacara. La mayoría de estas locaciones 
tiene como factor común estar ubicada a media 
altitud y relativamente próxima a la divisoria 
hidrográfica. Debe tenerse en cuenta que estas 
tendencias no son estadísticamente significa- 
tivas (cuadro 2); sin embargo, sugieren que el 
aporte pluviométrico en la zona montañosa es 
ligeramente más susceptible a presentar va- 
riaciones estacionales que su contraparte baja. 
Este sutil rasgo sugiere que el lago podría ser la 
principal fuente de humedad para la formación 
de lluvias en la parte baja de la cuenca, en parti- 
cular durante la temporada seca. Es importante 
destacar que es poco probable que los alisios 
que llegan a la cuenca del lago de Valencia 
estén cargados de abundante humedad, pues 
deben cruzar de manera previa las serranías de 
la costa (al norte) y del interior (al este). Ambas 
serranías bloquean el avance de los alisios (en 
especial, la serranía de la costa) y favorecen la 
ocurrencia de lluvias fuera del parteaguas de 
la cuenca (efecto Barlovento-Sotavento; deta- 
llado en Pulwarty et al., 1992). No obstante, los 
sistemas de tormentas de gran escala pueden 
ingresar a la cuenca, pues los niveles de bases 
de las nubes que lo conforman estarían sobre 
las cumbres más altas. 


Asociación entre la ocurrencia de extremos 
pluviométricos estacionales y las fases de 
ENOA en el periodo 1966-1992 


En la temporada seca del periodo 1966-1992 
(enero-abril), alrededor de 22% de los meses 
coincidió con un episodio La Niña; 19%, con El 
Niño, y 60% con la fase neutra. Al considerar los 
meses de la temporada húmeda (mayo-octubre), 
de manera curiosa esta categorización porcen- 
tual es idéntica. Por tanto, la incidencia de la 
fase fría de ENOA (La Niña) fue dominante en 
ambas temporadas. 

El cuadro 3 muestra que durante la tempo- 
rada seca a lo largo del periodo 1966-1992, El 
Niño coincide con la ocurrencia de un mes ES, 
EH y NE en la estación 489: 35, 15 y 50% de las 
veces, respectivamente, lo cual sugiere que los 
meses ES pueden atribuirse, de modo parcial, 
a la fase cálida de ENOA. De forma similar, en 
la estación 491 se encontró una asociación esta- 
dísticamente significativa entre las fases ENOA 
y la precipitación estacional (p < 0.05, basado 
en prueba x?). Sin embargo, en este caso, la fre- 
cuencia de ocurrencia de un mes ES o EH es más 
alta durante la fase neutra de ENOA. De hecho, 
un mes ES o EH muestra una probabilidad de 
ocurrencia de 20 y 17%, respectivamente, en au- 
sencia de El Niño o La Niña, sugiriendo que uno 
O varios factores no vinculados con ENOA están 
relacionados con estas condiciones climáticas a 
escala local. Al considerar la temporada húmeda 
(mayo-octubre), sólo resulta estadísticamente 
significativa la asociación ENOA-precipitación 
en la estación 1494 (figura 2). A diferencia de las 
estaciones anteriores, la respuesta en la estación 
1494 es dual, esto es, El Niño se vincula con 20% 
de los meses ES, mientras que La Niña coincidió 
con 25% de los meses EH. En general, tales re- 
sultados evidencian que la ocurrencia de un mes 
ES o EH a escala local durante las temporadas 
seca y húmeda no puede atribuirse sólo a las 
fases activas de ENOA, sino que otros factores 
juegan un papel más importante (probablemen- 
te factores climáticos locales). 

Durante la temporada seca (enero-abril) del 
periodo 1966-1992, se observó que, en prome- 
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dio, un mes ES, EH o NE es común en un 57, 
31 y 63%, respectivamente, entre las estaciones 
analizadas cuando se da un episodio El Niño; 32, 
18 y 73%, respectivamente, al suceder un episo- 
dio La Niña; y 39, 34 y 63%, respectivamente, en 
la fase neutra de ENOA (figura 3). A lo largo de 
la temporada húmeda (mayo-octubre), en pro- 
medio, un mes ES, EH o NE es común en 27, 21 
y 78%, respectivamente, del total de estaciones, 
al coincidir con un episodio El Niño; 33, 26 y 
76%, respectivamente, al suceder un episodio 
La Niña; y 24, 27 y 79%, respectivamente, en la 
fase neutra de ENOA (figura 3). Estos valores 
porcentuales revelan que los eventos El Niño 
durante la temporada seca se asocian con condi- 
ciones climáticas secas moderadamente extensas 
dentro de la cuenca (en promedio, 57%), mien- 
tras que la extensión espacial de los meses ES o 
EH durante un año La Niña es marcadamente 
menor a la observada durante los años El Niño. 
Al considerar la temporada húmeda, la inci- 
dencia espacial de El Niño es apenas superior 
a la observada durante la fase neutra (27 contra 
24%); mientras que La Niña, de forma curiosa, 
se asocia con meses ES y más extenso que los 
notados cuando coincide con El Niño (ES según 
cobertura: 27% para El Niño; 33% para La Niña; 
y 24% neutro). En la temporada seca del periodo 
1966-1992 (enero-abril) destacan dos eventos se- 
cos que cubrieron más de 80% de las estaciones 
en los años 1979 y 1987, vinculados con las fases 
neutra y El Niño de ENOA, respectivamente 
(figura 3a). Al considerar la temporada húmeda 
(mayo-octubre), sobresale la sequía del año 1973 
(figura 3b), la cual causó meses ES en más de 
60% de las estaciones y coincidió con la fase 
neutra de ENOA. 

El contraste de la figura 1 y el cuadro 3 
muestra que durante la temporada seca y años 
El Niño, los ríos Noguera y Las Minas (tributa- 
rios del río Giigie y la Laguna de Taguayguay) 
tienen una mayor probabilidad de presentar 
meses ES. La persistencia temporal de estas 
condiciones climáticas puede causar que el 
caudal de estiaje del río Gúigúe y el nivel de la 
laguna de Taguayguay disminuyan al tiempo de 
aumentar el riesgo de incendios forestales por 


la alta sequedad del suelo. En la temporada hú- 
meda de los años El Niño, los meses ES son más 
probables en la cuenca alta del río Las Minas y 
la cuenca baja del río El Limón; mientras que los 
meses EH son ligeramente más probables en los 
años La Niña sólo en Zuata y sus alrededores. 

Para el periodo analizado, los resultados des- 
critos evidencian que el fenómeno ENOA puede 
explicar de forma parcial la ocurrencia de meses 
ES o EH a escala local durante las temporadas 
húmedas y secas en la cuenca del lago de Valen- 
cia. El Niño presenta una mayor incidencia sobre 
los extremos pluviométricos que La Niña a esca- 
la local y suele estar relacionado con condiciones 
anómalamente secas con moderada cobertura 
espacial. La moderada a escasa conexión entre 
las fases activas de ENOA identificada a escala 
local sugiere que otros factores no evaluados 
incidieron de forma determinante sobre ciertos 
eventos pluviométricos extremos (por ejemplo, 
la sequía generalizada de 1979 se dio durante la 
fase neutra de ENSO). 

Cabe destacar que la asociación entre las 
fases ENOA y las lluvias en el territorio vene- 
zolano muestra en general una alta variabilidad 
espacial y temporal (Cárdenas et al., 2002). Esta 
característica ha sido atribuida a las complejas 
interacciones entre ENOA y ciertos modos de 
gran escala en el océano Atlántico y la varia- 
bilidad de la circulación atmosférica tropical 
(Martelo, 2003b). En todo caso, los resultados 
encontrados son coherentes con los hallazgos 
previamente reportados por Pulwarty et al. 
(1992), Martelo (2003b), y Paredes, Millano y 
Guevara (2008), entre otros, quienes señalan 
que los años El Niño favorecen la ocurrencia 
de meses ES; mientras que los años La Niña, 
meses EH. 


Conclusiones 


Durante los periodos 1934-2005 y 1966-1992 
no se evidenció un aumento o disminución de 
los acumulados pluviométricos estacionales 
en la cuenca del lago de Valencia. Aunque se 
identificaron tendencias locales de diferentes 
grados de magnitud, no son estadísticamente 
significativas. 
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A lo largo de las temporadas de seca y 
húmeda de los años El Niño o La Niña, sólo 
algunas locaciones muestran una probabilidad 
de ocurrencia de meses pluviométricos extremos 
(EH o ES) ligeramente mayor a la observada en 
años neutro. En general, las sequías vinculadas 
con El Niño suelen ser más extensas que las 
condiciones climáticas más húmedas que lo 
normal vinculadas con La Niña. En todo caso, 
una alta proporción de las sequías o temporales 
persistentes en la cuenca no están asociadas con 
las fases activas de ENSO, sino posiblemente con 
otros factores climáticos locales no identificados. 
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Anexos 


Prueba de autocorrelación serial de Durbin- 
Watson 


Al aplicar un análisis de regresión simple o múl- 
tiple entre dos o un grupo de variables se puede 
detectar la ocurrencia de una correlación seriada 
para un desfase unitario en los residuos deriva- 
dos del modelo ajustado. Suponga que la serie 
anual de precipitación en una estación k tiene 
L años de longitud, entonces sus observaciones 
forman un vector a,...a,, donde ¡= 1... L. Si las 
observaciones se ordenan en forma ascendente 
resulta p (x, y) = (1, a,), (2, a,)... (L, a,). Al aplicar 
un análisis de regresión lineal a la matriz p se 
obtiene una ecuación lineal ajustada, cuyos resi- 
duos forman un vector e (e,, e,, €,...e, ). Entonces, 
el estadístico de prueba Durbin-Watson (d) está 
dado por: 


a (1) 


El estadístico d es contrastado contra dos 
umbrales de. y Du anexado en Fox (2008). Estos 
umbrales dependen del nivel de significación 
considerado en la prueba (a), L y el número 
de predictores de la ecuación de regresión. Las 
reglas de decisión son: 


Si d < d, , los términos de error están auto- 
correlacionados de forma positiva. 

Si d > d,,, los términos de error no están 
autocorrelacionados de forma positiva. 

Si d, , <d<d,,,, la prueba no es concluyente. 

70 ¡0 

Si (4 —- d) < d, , los términos de error están 
autocorrelacionados de modo negativo. 

Si (4- d) > d,,, los términos de error no están 
autocorrelacionados de modo negativo. 

Si de < (4-d)< d,;, la prueba no es con- 
cluyente. 


La detección de una autocorrelación serial es- 
tadísticamente significativa en la serie anual de 
precipitación puede ser atribuida a la ocurrencia 
de errores en la medición de la precipitación. 


Prueba de homogeneidad de Easterling- 
Peterson 


Esta técnica permite detectar variaciones en la 
serie de precipitación anual inducidas por facto- 
res no climáticos. Suponga que la serie anual de 
precipitación en una estación k tiene L años de 
longitud de registros. En principio se selecciona 
una estación de referencia cercana a la estación 
evaluada con igual longitud de registro (al 
menos 30 años en común). Después se obtiene 
un vector diferencia, restando las observaciones 
de ambas estaciones año con año desde t= 1... 
L. Luego, dicha serie se divide en dos subseries 
S, y S, de longitud L, y L,, donde n =2... L; por 
tanto, L, =L—n, L,=L- L,, resultando S, =4,, 
Ayo A Y SR y A qua y=* A: Para identificar si 
el año n es no homogéneo, se aplica una prueba 
t-Student de diferencia de medias a las observa- 
ciones contenidas en cada subserie a un nivel de 
confianza aL. Si esta prueba resulta significativa, 
se concluye que la serie no es homogénea en el 
año 1. Una serie puede tener más de un año no 
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homogéneo, lo cual podría ser causado por un 
factor no climático. 


Prueba de tendencia de Mann-Kendall 


Esta prueba permite detectar una tendencia 
monótona en una serie de precipitación anual. 
Suponga que la serie anual de precipitación en 
una estación k tiene L años de longitud, enton- 
ces sus Observaciones forman un vector 4,... 
a, donde i = 1... L (note que L > 30 años). Las 
observaciones se analizan pareadas, como sigue: 
p (x, y) = (1, a,), Q, a,)... (L, a,). En principio, se 
calcula el parámetro S de Kendall: 

Para a, > si a,> a, score =+1, si a,< a, score=-1, 
si a,= 4, score = 0; si a,¿> a, score = +1, si a,< a, 
score =-1, si 4,= 4, score = O;... si 4, > 4, Score = +1, 
si a, < a, score =-1, si a, = a, score =0. Este proce- 
dimiento se repite desde a, hasta a,. Note que en 
cada contraste, 4,> 4,04,< a, siempre ¡> 1. 

Luego se totalizan los scores. P representa el 
total de scores positivos y M el total de los scores 
negativos. Hecho esto, S =P — M. El número to- 
tal de comparaciones está dado por L (L-1)/2. 
Por tanto, el 1 de Kendall es: 


S 
“ae? e 


Para estimar la significancia estadística del 
1 de Kendall, se consulta el anexo B de Helsel y 
Hirsch (1993), empleando los parámetros L y S 
como variables de entrada. 
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